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Abstract 



The classical approach to diffusion processes is based on Fick's law that the flux is pro- 
0h I portional to the concentration gradient. Various phenomena occurring during propagation of 

penetrating liquids in polymers show that this type of diffusion exhibits anomalous behavior 
and contradicts the just mentioned law. However, they can be explained in the framework 
of non-Fickian diffusion theories based on viscoelasticity of polymers. Initial-boundary value 
problems for viscoelastic diffusion equations have been studied by several authors. Most of the 
studies are devoted to the Dirichlet BVP (the concentration is given on the boundary of the 
domain). In this chapter we study the second BVP, i.e. when the normal component of the 
J> ' concentration flux is prescribed on the boundary, which is more realistic in many physical sit- 

uations. We establish existence of weak solutions to this problem. We suggest some conditions 
£T) '. on the coefficients and boundary data under which all the solutions tend to the homogeneous 

<0 state as time goes to infinity. 

t-H ■ 

g ; 1 Introduction 

The continuity equation for diffusion 

I--*,, (i.i) 

states that variations of the concentration u(t, x) at any spatial point x and moment of time t can 
only be caused by inflow and outflow of a penetrant into and out of that area. Here J = J(t, x) is 
the concentration flux vector. 

The classical diffusion theory is based on Fick's law (the flux is proportional to the concentration 
gradient with negative proportionality factor —D). The continuity equation and Fick's law yield the 
classical diffusion equation 

fin 

— = div{D{u)Vu), (1.2) 

which becomes the heat equation 

*= D A« (1.3) 

for constant diffusion coefficient D. 

The concentration behaviour for diffusion of penetrant liquids in polymers exhibits such phe- 
nomena as case II diffusion, sorption overshoot, literal skinning, trapping skinning and desorption 
overshoot, which contradict Fick's law, see e.g. [SJ El EJ EH EDI ED E31 ES]- There is a number 
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of approaches which explain these non-Fickian properties of polymeric diffusion. They are usually 
based on taking into account the viscoelastic nature of polymers (cf. [15] and references therein) and 
on the possibility of glass-rubber phase transition (see e.g. [26] with some review). We are going to 
study the model which is due to Cohen et al. (see [U O [11] ; related models or particular cases of 
this one were suggested by other authors, see e.g. [H It consists in combining the continuity 

equation with the following system 

J(t,x) = -DVu- EX7a + Mu, (1.4) 

— + (3a = [iu + vu . (1.5) 

Generally speaking, the coefficients (3, D, E, M, \i and v may depend u, or even on t, x and a. 
Let us briefly discuss the meaning and typical behaviour of the coefficients. The scalar function 
(3 is the inverse of the relaxation time. A characteristic form of (3 is [5] 

(3(u) = k(3 R + (3 G ) + ~(f3 R - Pa) tanh(^-^) (1.6) 
ll o 

where (3r, (3g,5,u rg are positive constants, (3r > (3g- The polymer-penetrant systems modeled with 
the help of (1.6) can be in two phases: glassy and rubbery. The glassy state corresponds to the areas 
of low concentration. Here the polymer network is severely entangled, and the relaxation time is 
high, so its inverse is low. Moreover, it is close to a certain value (3g- In the high concentration areas 
the system is in the rubbery state: the network disentangles, so the relaxation time is small, and its 
inverse is close to (3r > (3g- The glass-rubber phase transition occurs near a certain concentration 
urg, and the value of 5 determines the length of the transition segment. The coefficients D and E are 
non-negative scalars (more generally, they are positive-definite tensors) called the diffusion and stress- 
diffusion coefficients, respectively. As the concentration increases, the polymer network disentangles, 
so the diffusivity also increases. Thus, D should be an increasing function of concentration: in 
particular, D can depend on u in a way similar to (1.6) [10]. E is sometimes considered to be a 
constant, see [I] for some justification, but numerical simulations [6] have shown that, if E(0) ^ 0, 
then the concentration u may become negative, which is physically meaningless. Conversely, it can 
analytically be proved that, if E(0) = 0, then the concentration u remains non-negative provided 
it is non- negative at the initial moment of time [3]. In [23], we make related observations showing 
expediency of the condition E(l) = 0, which can maintain the concentration u of less than or equal 
to 100%. Thus, a modeling example is 

E ( u > = n ?\2' ( L7 ) 

where a\ and a 2 are positive constants, a 2 is small. The functions \i and v should be non-negative 
and bounded [3], M is the convection velocity vector, assumptions on it will be given below. 

Let n(x) be the outward normal vector at the point x of the boundary dfl of a domain Q C M" - 
Q. Then system (1.1), (1.4), (1.5) may be completed with such boundary conditions as 

u(t,x) = <p(t,x), xedVt (1.8) 

(the concentration on the boundary is prescribed) and 

n 

— Jj(t, x)nj(x) — cp(t, x), x G dQ (1.9) 
i=i 



lr The most important particular cases are n = 2 (diffusion in polymer films) and n = 3. 
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(the influx^ of the penetrant through the boundary is known). 

The initial-boundary value problems for system (1.1), (1.4), (1.5) possess maximal (not global in 
time) solutions for a more general boundary condition, which includes (1.8) and (1.9), see [3]. The 
global (in time) existence results are known for the Dirichlet condition (1.8). A theorem on global 
solvability is presented in [2j for / = uu, M = and D = E being a constant scalar. It is formulated 
for the one-dimensional case (0 < x < 1), but the technique used there seems to be applicable for 
x £ Q, where Q C R n is a bounded domain with a smooth boundary. Another global existence 
result is given in [13]. They assume the stress-diffusion coefficient E to be a non-constant increasing 
function of concentration, E(0) = (so E(l) = is not allowed). However, it is required that the 
initial and boundary data for the concentration are bounded from below by a positive constant, so 
the solution is always strictly positive, and this approach does not permit to consider dry regions 
in a polymer. Paper [13] is mainly concerned with the one-dimensional case, but also they suggest 
a brief plan how to generalize the result for the multidimensional situation. Global existence of 
dissipative (ultra weak) solutions for constant scalar D and E and M = is shown in [22] for 
Q = R™ (again, the ideas used there seem to be suitable for Q C R"). In [23], global (in time) weak 
solutions on a bounded domain Q C R n are constructed, under rather general assumptions on the 
coefficients. Further investigation of the weak solutions of the Dirichlet problem is carried out in |24j : 
it is proved that, for any sufficiently short time segment and any stress prescribed at the beginning 
of this segment, there exists a weak solution such that the concentrations at the beginning and at 
the end of the segment are the same, and, under an additional assumption on coefficients, existence 
of time-periodic weak solutions (without any restrictions of the period length) is shown. Paper [23] 
considers long-time behaviour issues for this problem: provided D and E are constant scalars and 
M = 0, the solutions generate a dissipative semiflow, and there exist a minimal trajectory attractor 
and a global attractor. 

In this chapter, we construct (global in time) weak solutions for problem (1.1), (1.4), (1.5), (1.9) on 
a bounded domain Q C R™ for given initial concentration and stress. The coefficients may depend on 
t, x, u and a. In addition, we suggest some conditions on the coefficients and boundary data under 
which all the solutions tend to the homogeneous state u = const as time goes to infinity. The chapter 
is organized in the following way. In Section 2, we introduce the required notations. In Section 3, 
we give a weak formulation of the initial-boundary value problem and state the result on existence 
of weak solutions (Theorem 3.1), which is proved in Section 4. In Section 5, we touch the long-time 
behaviour. 

2 Notation 

We use the standard notations L P (Q), W™(Q), H m (Q) = W™(fl) for Lebesgue and Sobolev spaces 
of functions defined on a bounded open set (domain) Q C R™, n G N. 

The scalar product and the Euclidean norm in L2{Q) k = L2(fi,R fc ) are denoted by (u,v) and 
||w||, respectively (k is equal to 1 or n). In H m (Q), m e N, we use the scalar product (u,v) m = 
{D a u, D a v) and the corresponding Euclidean norm ||u|| m . 

|a|<m 

The space of linear continuous functionals on H m (Q) (the dual space) is denoted by H^ m (Q). 
The value of a functional from H^ m (Q) on an element from H m (Q) is denoted by (•, •) (the "bra-ket" 
notation). Similarly, the dual space of W™(Q) is denoted by W~™(£1), ~ + ~ = 1, 1 < p < oo, with 
the corresponding use of the "bra-ket" notation. 

Sometimes we shall write simply L p , H m for L p (Q) k , H m (fl) k etc., k = 1, n. 

2 It can be negative. 
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Let us introduce some basic operators. The operator div^ : L q (Q) — > W ^(Q) is determined by 
the formula 

f 11 

(div N v, 4>) = - v(x)V(j)(x) dx, <\> G W*(tt), - + - = 1. (2.1) 

J p q 

n 

The isomorphic operators A : if 1 (ft) -> H^^Q) and A 2 : 22 2 (ft) -> 22~~ 2 (ft) are given by the 
expressions 

(Au, 0) = (f , 0)i, G 22 1 (ft), (2.2) 

(^ 2 u,0) = (t; j 0) 2 , </>G 22 2 (ft). (2.3) 

Note that L 2 (ft) C 22^ (ft) C 22j^ 2 (ft) with natural imbedding operators, and then = v — div^^v, 
v G H 1 . 

Set X N = X N (Q) = A _1 (if 1 (f2)). The scalar product and norm in X N are (u, v )x = (Au,Av)i, 
\\u\\ x = ||Au||i. The duality between 22^* (ft) and X N (fl) is given by the formula 

(u,v) 1 = (u,Av) , u e H^\v e X N . (2.4) 

Note that (u,v) 1 = (u,v)i for u G H l ,v G Xjv- The elements of X^ are solutions of the Neumann 
problem 

v- Av = ue H l (tt), (2.5) 
dv 

— (x) = 0,xedn. (2.6) 
on 

Thus, 

x N (n) c 22 2 (ft) c wl /n _ 2 {tt) (2.7) 

by Sobolev theorem (for sufficiently regular ft). 

The symbols C(J; E), C W {J\ E), L 2 (J; E) etc. denote the spaces of continuous, weakly contin- 
uous, quadratically integrable etc. functions on an interval J C R with values in a Banach space E. 
We recall that a function u : J — > E is weakly continuous if for any linear continuous functional g 
on E the function g(u(-)) : J — > R is continuous. 

If 22 is a function space (L 2 (ft), 22" m (ft) etc.), then we identify the elements of C(J; E), L 2 {J\ E) 
etc. with scalar functions defined on J x ft according to the formula 

u(t)(x) = u(t,x), t G JT, x G ft. 

We shall also use the function spaces (T is a positive number): 

W N = W N (Q,T) = {re L 2 (0,T;H\Q)), r' G L 2 (0,T; H N \Q))} 

\\ T \\w N = \\r\\L 2 (0,T;H^(n)) + lk'||L 2 (0,T;H^ 1 ( n ))' 

Wi = = {r G L 2 (0,T;X w (fi)), r' G L 2 (0,T; H^(0))} 

ll T lki = ll r IU 2 (0,T;X i v(n)) + \W\\l 2 (0,T;H^ 1 (U)^ 

W 2 = W 2 (Q,T) = {re L 2 (0,T;22 2 (ft)), r' G L 2 (0, T; 22~ 2 (ft))} 

\\ T \\W 2 = ||T||jr ia ( ,T;fla(n)) + \\ T '\\L 2 {0,T;H- 2 (n))- 

Lemma III. 1.2 from [19] implies continuous embeddings Wn,W 2 C C([0, T]; L 2 (ft)), W± C 
CdO^H 1 ^)) (see also [L2]). 
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We use the notation | • | for the absolute value of a number, for the Euclidean norm in R" - , and 
in the following case. 

Denote by W nxn the space of matrices of the order n x n with the norm 

|Q| = max \Q£,\. 

£6RMI|=1 

Let R" x ™ Cl" x " be the set of such matrices Q that 

for some d(Q) > and all f G R n . 

The symbol C will stand for a generic positive constant that can take different values in different 
places. 

3 Weak formulation of the problem 

We consider a polymer filling a sufficiently regulajf] bounded domain O C I", n G N. We study the 
diffusion of a penetrant in this polymer which is described^ by the following initial-boundary value 
problem^]: 

du 

— = mv[D (t, x, u, a) Vw 

+E (t,x,u,a)Va - M Q (t,x,u,a)u], (t,i) 6 [0,T] xO, (3.1) 
da du 

— + (3o(t, x, u, a)a = u (u)u + v (u) — , (t, x) e [0, T] x O, (3.2) 



du „ „ s da 



+ E (t,x,u,a)ij- 



, X, It, I7J jj | ~ yjj 



-M (t,x, u,a) l u n.i(x) = <p{t,x), (t,x) G [0,T] x dQ, (3.3) 

w(0, ar) = uo(x), a(0,x) = a (x), x E Q. (3.4) 

Here u = u(t,x) : [0,T] x — > R is the unknown concentration of the penetrant (at the spatial 
point x at the moment of time t), a = a(t,x) : [0, T] x Q — > R is the unknown stress, w = Mo(x), 
°"o = cr o( a; ) : ^ ~~ K are given initial data, (p : [0,T] x <9f2 — > R is the influx of the liquid through 
the boundary, u ,u : R -> R, D , £ : K n+3 = RxR n xRxR^ R™ xri , (3 : R n+3 -> R, 
M Q : R n+3 — ► R n are given functions, u (-) > 0. 

Before giving a definition of a weak solution to this problem, it is convenient to make a change 
of variables. Denote 

u(t,x) 

q(t, x) = a(t, x) - I u (y) dy, 



o 



3 Say, it is locally located on one side of its C 2 -smooth boundary. 

4 This problem makes sense for diffusion in polymers provided < u < 1, i.e. the concentration is not less than 0% 
and does not exceed 100%. The assumptions on coefficients which guarantee this condition provided it is fulfilled at 
the initial moment of time are discussed in (24j for the Dirichlet problem, mainly based on the results of [3]; similar 
arguments are applicable in the Neumann case (1.9) with ip = 0. However, this question is still not completely studied. 
So we consider here the general setting (3.1) - (3.4). 

5 System (3.1), (3.2), (3.3) is obtained from (1.1), (1.4), (1.5), (1.9). For technical purposes, we assign subscript zero 
to the coefficients. 
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uo(x) 

q (x) = a (x) - J u (y) dy, 
o 

D(t,x,u,q) = 



D lt,x,u,q + J u (y)dy J + vq{u)E q j t,x,u,q + J v (y)dy J G 



artxn 







it, <r) = £ j t, x, u, q + J u (y)dy J , 



o 

u 



f(t, x, u, q) = -uM lt,x,u,q + J v (y)dy J , 



o 

u 



/?i(t, x, u, = -fa lt,x,u,q + J v (y)dy J , 



j(t,x,u,q) = (io(u) - 



u \ u 

po ( t, x, u, q + J v (y)dy J v (y) dy 
o / o 



u 

Note that, if u vanishes, then, by continuity, we consider the last term to become 

-P o (t,x,0,q)v o (0). 
Then we can rewrite (3.1) - (3.4) in the following form: 

du 

— = div[D(t, x, u, q)Vu + E(t, x, u, q)S7q + f(t, x, u, q)}, (3.5) 
dt 

dq 

— = (3i(t, x, u, q)q + 7(t, x, u, q)u, (3.6) 

Er . du dq 

\D(t, x, u, VijQ^ + E (t, x, u, VijQ^ 

i,j=l 3 3 

+f(t, x, u, q)^ Hi(x) = (p(t, x), (t, x) E [0, T] x dQ, (3.7) 

u\ t=0 = w , q\t=o = <S>- (3.8) 

Now, before describing our assumptions on the coefficients, let us calculate the gradient of the 
right member of (3.6): 

V(/?i(£, x, u, q)q) + V{~f(t, x, u, q)u) 

= Pi(t, x, u, q)Vq + ^-(t, x, u, q)q + ^-(t, x, u, q)qVu + ^-(t, x, u, q)qVq 
dx du dq 

d^j d^ d^y 

+j(t, x, u, q)X7u + — (t, x, u, q)u + — (t, x, u, q)uVu + — (t, x, u, q)uVq 

= (3(t, x, u, q)Vu + a(t, x, u, q)Vq + g(t, x, u, q), (3.9) 

where 

d(3i d-f 

/3{t, x, u, q) = -Q^{t, x, u, q)q + j(t, x, u, q) + —(t, x, u, q)u, (3.10) 
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(i,(t, x, u, q) = f3i(t, x, u, q) + -Q-{t, u, q)q + -Q-(t, x, u, q)u, (3.11) 

g(t, x, u, q) = -^(t, x, u, q)q + ^-(t, x, u, q)u. (3.12) 

We assume the following: 

i) D, E : M n + 3 -> R nxn ; f, g : M n+3 -> W 1 ; //, /3, 7, ft : M n+3 -> R. 

ii) Each of these eight functions (e.g. D(t,x,u,q)) is measurable in (t,x) for fixed (u,q). 

iii) Each of these functions is continuous in (u,q) for fixed (t,x). 

iv) These functions satisfy the estimates 



\D(t,x,u,q)\ < K D , 


(3.13) 


\E(t,x,u,q)\ < K E , 


(3.14) 


max(\f3(t,x,u,q)\, \j(t,x,u,q)\) < Kp, 


(3.15) 


raax(\fi(t,x,u,q)\, |/3i(t, x, u, q)\) < K^, 


(3.16) 


\f(t,x,u,q)\ < K f {\u\ + \q\) + f{t,x), 


(3.17) 


\g(t,x,u,q)\ < K g (\u\ + \q\)+g(t,x) 


(3.18) 



with some constants K D , . . . , K g and functional f,g G L 2 i oc (W n+1 ). 
v) 

(D{t,X,U,q)q-,ORn>d(q-,0 R n, (3.19) 

where d > is independent of (t, x, u, q) G M n+3 and £ G M n . 
vi) Relations (3.10) - (3.12) hold. 

It is easy to see that, if -Bo and (3q are taken in the forms (1.7) and (1.6), then (3.14) and (3.15) are 
violated. It turns out that such deficiencies can be corrected without loss of generality of the model 
(see j2H Section 3] for a detailed discussion, cf. also [231 125]). In brief, physically, the concentration 
u and the stress q are uniformly bounded, therefore the coefficients of systems (3.1)-(3.2) (and, 
consequently, of (3.5)-(3.6)) can be experimentally determined only for bounded u and q, whereas 
"at infinity" we can choose them at discretion. 

Let us now rewrite (3.5) and (3.7) in a weak form. Assuming u and q sufficiently regular, take 
the L 2 (n)-scalar product of the members of (3.5) with a test function <fi G H 1 ^), and integrate by 
parts in the right-hand side: 



+ ^ / [D(t,x,u,q)ij^ + E(t,x,u,q)ij^ + f(t,x,u,q)i ni(x)<f)(x)ds 



(if', 4>) = -(D(t, x, u, q)Vu + E(t, x, u, q)Vq + f(t, x, u, q), V</>) 

dq_ 

'dxj 

= -(D(t,x,u,q)Vu + E(t,x,u,q)Vq + f{t,x,u,q),V<f)) + J <pcf>ds. (3.20) 

an 

Denote by ip(t) the linear functional (p *—>■ J (p(t)<pds. We assume that (for a. a. t) this integral exists 

an 

and continuously depends on <p G H 1 ^), so ip(t) G H^ 1 ^); clearly, this is true e.g. if <p(t) G L 2 (d{l). 
Then we arrive at 

Ou 

— = div N [D(t, x, u, q)Vu + E(t, x, u, q)Vq + f(t, x, u, q)] + if), (3-21) 



'Clearly, the behaviour of these functions outside (0, T) x Q does not matter. 
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which should be understood as an equality of functionals from i/ 7V 1 (f2). Conversely, for each pair of 
sufficiently regular functions (u,q), (3.21) implies (3.5) and (3.7). 

Definition 3.1. A pair of functions (u,q) from the class 

u G W N (n,T),q G tf^T;// 1 ^)) (3.22) 

is a weak solution to problem (3.5)-(3.8) if it satisfies (3.8), equality (3.21) holds in the space i/ 7 ^ 1 (f2) 
a.e. on (0,T), and (3.6) holds a.e. in (0,T) x Q. 

Note that (3.8) makes sense due to the embeddings 

W N C C([0,T];L 2 (Q)), H\0, T; H\n)) C C([0, T\- H\Q)). 

Theorem 3.1. For every u G L 2 (Vi), % G H 1 ^) and ip G L 2 (0, T; H^ 1 ^)), there exists a weak 
solution to problem (3.5) - (3.8) in class (3.22). 



4 Proof of the existence result 

The proof of Theorem 3.1 is based on the study of the following auxiliary problem: 

dv 

— + sA 2 v = Xdiv N [D(t, x, v, r)Vv + E(t, x, v, r)Vr + f(t, x, v, r)] + A^, (4.1) 
dr 

— + sA 2 t = A[/?i(t, x, v, t)t + 7 (t, x, v, t)v], (4.2) 

v\ t=0 = u , (4.3) 

r\ t =o = (4.4) 

Here e > 0, A e [0, 1] are parameters. We are going to derive some a priori estimates for the weak 
solutions of this problem. Then we shall show its solvability via topological degree arguments (the 
presence of the parameter A is important at this stage). Finally, we shall put A = 1 and pass to the 
limit as e — > 0. 

Definition 4.1. Given u G L 2 (Q), q G if 1 (fi), a pair of functions (v,t) from the class 

v G W 2 (n,T),r G W^T) (4.5) 

is a weak solution of problem (4.1)-(4.4) if equality (4.1) holds in the space H^ 2 {VL) a.e. on (0, T), 
(4.2) holds in the space H N l {tt) a.e. on (0,T), (4.3) holds in L 2 (Q), and (4.4) holds in H\Q). 

The last two conditions make sense due to the embeddings 

W 1 C C([0,T];H\Q)), W 2 C C([0, T}; L 2 (Q)). 

Lemma 4.1. Let (v,t) be a weak solution to problem (4-l)-(4-4)- Then the following a priori 
estimate holds: 

1 1 1 1 2 1 1 1 1 2 

£ \\ V \\L 2 (0,T;m{n)) + £ \\ T \\l 2 (0,T;X n ) + 

ll v llioo(o,T;L 2 (n)) + M\ v \\L 2 (o,T;m(n)) + IMILoCo.tw 1 ^)) - C ( 4 - 6 ) 
where C is independent of A and e. 
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Proof. Take the "bra-ket" of the terms of (4.2) (as elements of ^(fi)) and Ar(t) G H 1 ^) at 
a.a. t G [0,T]: 

(r', At) + (eA 2 t, At) 
= X (J3i(t, x, v, t)t + j(t, x, v, t)v, At) . (4.7) 
Note that we can use parentheses instead of brackets in the right-hand side due to the equality 

(wi,w 2 ) = (wi,w 2 ), Wi G L 2 ,w 2 G H l . 

But 

(r',AT) = (T',T) 1 = 1 -j t \\T\\l (4.8) 
(e.g. by PSI Lemma III. 1.2]). Thus, 

= X (Pi(t, x, v, t)t + j(t, x, v, t)v, t) + A (V[(3i(t, x, v, r)r + j(t, x, v, r)v], Vr) 
= A x, r)r, r) + A (7(4, x, v, t)v, t) 
+A x, v, t)Vv + x, v, t) Vr + x, t>, r), Vr) . (4.9) 
Denote v(t) = e~ kt v(t), f(t) = e _fc V(t), where k > will be defined below. Then 



\(/3i(t, x, e kt v(t), e kt f(t))fe kt , r(t)e fct ) + \U(t, x, e kt v(t), e kt f (t))ve kt , f(t)e k 
+A x, e fc *t;(t), e fct f (t))Vve H + /i(t, x, e kt v{t),e kt f(t))Vfe kt 

+g(t, x, e kt v(t),e kt f(t)),Vf(t)e kt ) . (4.10) 



Denote now 



f3 k (t, x, v(t),T(t)) = /3(t, x, e kt v(t), e kt f (t)), 
Hk(t, x, v(t), f(t)) = (i(t, x, e kt v(t), e kt f{t)), 
g k {t, x, v(t), f (t)) = e- fc ^(t, x, e fct t;(t), e fci f (t)), 
ft fc (t, x, i;(t), f (£)) = A(t, x, e fci t;(t), e fct f (t)), 
7fc (t, x, i;(t), f (£)) = 7 (t, x, e fct t;(t), e fci f (*)). 



Thus, 



n^\\ f \\l + k \\ f \\l + £ ( f > f )x 



l d 

2~dt 

A | lk (t, x, f (t))f (t), f (t)J + X(j k (t, x, v(t), f(t))v(t). f( I ) 
+A (j3 k (t, x, v(t),f(t))Vv + /i k (t, x, v(t), f(t))Vf 

+^(t,x^(t),f(t)),Vf(t)Y (4.11) 
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Integration from to t G [0,T] yields 

t t 
+ * / \\f(s)\\lds + e J \\f(s)\\ 2 x ds 



t 

= ^lko||? + A J ((3 k (s,x,v{s),f{s))Vv(s),VT(s)) + (-f k (s,x,v(s),f(s))v(s),T(sf) 



+ (jj, k (s,x,v(s),f(s))Vf(s), Vt(s)J + ^i fc (s,x,u(s),f(s))f(s),f(s) 



+ ^ fc (s,a;,u(s),f(s)), Vf(s)) ds. (4.12) 

Applying the Cauchy-Buniakowski inequality, Cauchy's inequality ab < ca 2 + j-b 2 (3.15) and (3.16) 
we obtain 



l\\m\i + k J ||f( a )||?d a + e J \\f(s)\\ 2 x ds 



o o 
t t 



1 A./v CP f 

< 2NI? + -T ./ H«(*)ll? ds + A y \\T(s)\\lds + \K,j \\f(s)\\lds 



t t 
+^ | M S ,-^( S ,.),f( S ,-))H 2 ^ + A y ||f( S )||?da. (4.13) 



Note that 



||^(s,.,t;( S ,.),r(3,.))|rds< / -)l + \f(s, -)|] + -)lr ^ 



o 

i t 

■)ll 2 ds 



<3^y \\v(s,-)\\ 2 ds + 3K 2 J \\f(s,.)\\ 2 ds + 3 J \\g(s, 



t t 



-2 




< 3K 2 I \\v(s 



J \\v{s)\\\ds + 3K 2 g J ||f(a)||?d S + 3||flf||i a((0i r )xn) . 



Hence, 



l -\\f{t)\\\ + {k-2-K,-\K 2 g ) J \\T(s)\\lds + e J \\f(s)\\ 2 x ds 



<^IMI? + A(^ + ^ 2 ) J ||^)|| 2 ^ + ^||5llL( (0) T)xn)- (4-14) 



o 

t 



Take k > 4 + 2K^ + \K 2 g . 



In particular, (4.14) implies 

t t 
J \\f(s)\\lds<^(l + X J Ms)\\lds). (4.15) 



o 
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Now, take the "bra-ket" of (4.1) (as elements of i^ 2 (fi)) and v(t) G H 2 {9) at a.a. t G [0,T]: 

(?/, v) + (eA 2 v, v) 

= A (div N [D(t, x, v, t)Vv + E(t, x, v, r) Vr + f(t, x, v, r)\ +ip,v). (4.16) 
Again, by [191 Lemma III. 1.2], 



Thus, 



\v\\ 2 + e(v, v) 2 



1 d 
2dt' 

= -\(D(t, x, v, t)Vv + E(t, x, w, r)Vr + /(£, x, v, t),Vv) + A v) . (4.18) 

Then 

= -A(z>(t, x, e fci ^(t), e fct f (t))Vt;e fci + x, e fci ^(t), e fct f (t))Vfe fct 

+/(*, x, e kt v{t), e ht f{t)),Vv{t)e ht \ + A <^(t), v{t)e kt ) . (4.19) 



Denote now 



Thus, 



D k {t, x, v(t),f(t)) = D{t, x, e kt v(t),e kt f(t)), 
E k (t, x, v(t),f(t)) = E(t, x, e kt v{t), e kt f (t)), 
f k (t, x, v(t),f(t)) = e- kt f{t, x, e kt v(t), e kt f(t)). 



lj t \M 2 + k\M 2 + <v,v) 2 



Therefore 



= -\\D k (t, x, v(t),f(t))Vv + E k (t, x, v(t), f(i))Vf 

+/ fc (t, x, «(t),r(t)), Vv(*)) + e" fct A (^(t), «(*)) . (4.20) 

t t 
-\\v(t)\\ 2 + k [ \\v(s)\\ 2 ds + e f \\v(s)\\ 2 2 ds 



t 

2 



A / (D k (s,x,v(s),T(s))Vv(s) + E k (s,x,v(s),T(s))VT(s) 







+f k (s, x, t5(s), f (a)), Vu(s)) - e" fcs (V>00, *>(*)> ds. (4.21) 
Using Cauchy's inequality, (3.14) and (3.19), we get 



-\\v(t)\\ 2 + k J \\v(s)\\ 2 ds + e J \\v(s)\\ 2 ds + Xd J (Vv(s),Vv(s))ds 



t t 
<^Kf + ^ / \\Vf(s)\\ 2 ds + ^ J \\Vv(s)\\ 2 ds 
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{ { 

+2 I \\f k (s,;v(s,-),f(s,-))\\ 2 ds + ^ J \\Vv(s)\\ 2 ds. 



t t 
A f ,, , „ „ , |2 Ao /" ... . ||2 



+ ^ / ll#0ir-i^ + T / 11*0011! (4-22) 



As for g fc above, we have 



J \\fk{s,-,v(s,-),f( 



s,-)Wds 



t 



<3K 2 f J \\v(s)\\ 2 ds + 3Kj J ||f( S )||^ S + 3||/||| 2((0;T)xn) . 



Hence, from (4.22) and (4.15), 

2 t t t 

\m\ 2 + (k- 3 -^-±) J Ms)fds + e J \\v(s)\\lds + ^ J \\Vv(s)fds 



2 



< - 
~ 2 





t 

< y(l + A J \\v(s)\\lds)+C 



t t 

= ^(1 + A J \\v(s)\\ 2 ds + \ j \\Vv(s)\\ 2 ds) + C. (4.23) 



Take k> -^ + ^ + ± + f + l. Then (4.23) yields 

t t 
J \\v(s)rds+^ J \\Vv(s)fds<C 







(now C may depend on k), so 

•/ 11*0011? 



t 

A y ||u(s)||?ds < C. 
o 

Thus, the right-hand members of inequalities (4.14) and (4.23) are bounded, and we arrive at 



1 1 — 1 1 2 I ' — 1 1 2 

£ IMlL 2 (0,T;ff2(Q)) + £ II T \\l 2 (0,T;X n ) + 
ll^llioo(0,T;L 2 (n)) + ^ll^lli 2 (0,T;if 1 (f2)) + II f II LooiO^iH 1 (Q)) ^ @- (4.24) 

Since e kt < e kT for t E [0,T], this implies (4.6). □ 



12 



Lemma 4.2. Let (v,t) be a weak solution to problem (4-l)-(4-4)- Then there is the following bound 
of the time derivatives: 

W V '\\L 2 (0,T;H N \n)) + \\ T '\\L 2 {0,T;H-\n)) < C ( l + V^) (4.25) 

where C is independent of X and e. 

Proof. Really, since H^ l (VL) C H^ 2 (Vl) continuously, (4.1) and (4.6) imply 

W V '\\L 2 (0,T;H N 2 (n)) < £ \\ A 2 V \\ L 2(0,T;H N 2 (n)) + 

X\\div N [D(t, x, v, r)Vv + E(t, x, v, r) Vr + /(*, x, v, r)] || L2{0jT . H -^ n)) + A||^|| l^t-h-^q)) 

< VeV£\\v\\ L2 (o,T;H2(n)) + 
X\\D(t, x, v, t)Vv + E(t, x, v, t) Vr + /(*, x, v, r) |U 2 ( ,r ; L 2 (n)) + IMI ^(o^if- 1 ^)) 

< Cv^ + ^DA||v|| L2(0 ,T;ifi(n)) + ^E||T||L 2 (o,T;ifi(Q)) + ||/(*, Z, V, t) \\l 2 (0,T;L 2 (Q)) + C 

< C\fe + K D \fX\v\ L2 {§^-,m{VL)) + -^IMIl^tw 1 ^)) 

+-^/IMU 2 (0,T;L 2 (n)) + ^/||r||L 2 (0,T;L 2 (n)) + ||/|U 2 ((0,T)xfi) + C 

< C\fe + C[VX\\v\\ L2 (o :T]H i(n)) + ||T|| Loo ( 0) r;i?i(Q)) 

+ lblUoo(0,T;L 2 (n)) + IMUoc^TW 1 ^)) + 1] < C(l + V^)- 

Similarly, since L 2 (fi) C ^(fi) continuously, (4.2) and (4.6) yield 

ll r 'llL 2 (0,T;H^ 1 (n)) ^ £ ll^ 2r llL 2 (0,T;/i'^ 1 (n)) + 

\\\Pl(t, X, V, t)t + 7(t, X, v, r)v|| L2(0 ,T;L 2 (n)) 

< V / £V / i|klU 2 (0,T;X JV ) + 
^ll T IU2(0,T;L 2 (n)) + ^/3|k||L2(0,T;L 2 (n)) < C(l + . 

□ 

Lemma 4.3. Given u G L 2 (fi), ?o G -^ 1 (^); there exists a weak solution to problem (4-l)-(4-4) ^ n 
class (4-5). 

Proof. Let us introduce auxiliary operators by the following formulas: 

Q 1 :W 2 xW 1 ^L 2 (0,T;H N 2 (n)), 

Qx(v, t) = div N [D(-, v, t)Vv], 
Q 2 :W 2 xW 1 ^L 2 (0,T;H N 2 (n)), 

Q2(v, t) = div N [E(-, •, v, t)Vt], 
Q 3 :W 2 xW 1 ^L 2 (0,T;H N 2 (n)), 

Qs(v, t) = div N [f(-, v, r)] + ip, 
Q 4 :W 2 xW 1 ^L 2 (0,T;H N \Q)), 

Q 4 (v,T) = 7(-, ;V,T)v, 

Q 5 :W 2 xW 1 ^L 2 (0,T;H N 1 (n)), 
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Q 5 (v,t) = /3i(-, -,v,t)t, 



Q : W 2 x W x -> L 2 (0, T; x L 2 (0, T; ^(fi)) x L 2 (fl) x 

<5(u, r) = (-Qi(f , r) - Q 2 (u, r) - Q 3 {v, r), -Q 4 (u, r) - Qs{v, r), 0, 0), 

Ax : Wi -> L 2 (0, T; ff^(fi)) x ff 1 ^), 

ii(«) = (u' + eA 2 u,u\ t =o), 
A 2 :W 2 ^ L 2 (0,T;H N 2 (n)) x L 2 (Q), 
A 2 {u) = (u' + eA 2 u, u\ t =o), 
A:W 2 xW 1 ^ L 2 (0, T; i/^ 2 (fi)) x L 2 (0, T; x L 2 (Q) x F 1 ^), 

i(v,r) = {v' + eA 2 v,r' + eA 2 T,v\ t= o,T\ t=Q ). 
Then the weak statement of problem (4.1) - (4.4) is equivalent to the operator equation 

A(v, r) + XQ(v, r) = (0, 0, u , To)- (4.26) 

Let us briefly explain the idea of the proof. We are going to show that the operator A is invertible. 
This yields the solvability of equation (4.26) for A = 0. On the other hand, Q turns out to be a 
compact operator. Then we can rewrite (4.26) in a form suitable for application of the Leray-Schauder 
degree theory, which will imply the existence of solutions for all A G [0, 1]. 

We recall that a non-linear operator K : X\ — > X 2 (X\ and X 2 are Banach spaces) is called 
compact if it is continuous and the image of any bounded set in X\ is relatively compact in X 2 . 
In particular, if X\ is reflexive, and, for any sequence x m — > x* which converges in X\ in the weak 
sense, one has K(x m ) — > K(x*) strongly in X 2 , then K is compact (since any bounded subset of Xi 
is relatively compact in the weak topology). 

For some q > 2, the embeddings W\ C L,(0, T; W\(Sl)), W 2 C L q (0, T; ^(fi)) are compact. 
Really, we have W x C C([0, T]; iJ 1 (f2)) ; W 2 C C([0, T]; L 2 (0)) continuously. Note that (by the 
Rellich-Kondrashov theorem) H 2 C L 2 compactly. Furthermore, H 1 C -f/^y 1 compactly, so the adjoint 
embedding X N C i/ 1 is also compact. Then, by [TH Corollary 6], W x C L p (0, T; H x (££)), W 2 C 
L p (0, T; L 2 (Q)) compactly for every p < oo. Let Pl = ^ and Po = cf. (2.7). Then |- = | + i. 
For -u G Xat, we have 

Nl 2 ^ < C7(|HH pi + ||v M ||I w ). 

The second term is 

|||Vw| 2 || Lpi/2 < ||Vm||l 2 ||Vm|| Lpo < (7||«||i||m||x- 

The first term can be estimated similarly. If q\ > 2 is such that — = \ + - with some p large enough, 
then, by [TTJ Lemma 11], W\ C L qi (0,T; W x ) compactly. Now, for u G H 2 , we have 

ll«llw* < c(\\u\\i pi + \\vu\\i n ), 



and the second term is 

I ; ~ - I \ "11/ I \ "II; 



IVwl 2 !! 2 - < IIVull? JlVd' 2 



< C||m|| 2 ||m|| 2 < C||m||||m||^. 

The last inequality follows from [1], Theorem 4.17]. The first term can be estimated in the same 
way. If q 2 > 2 is such that ^ = | + - with some p large enough, then, by [T7J Lemma 11], 
W 2 C L q2 (0, T; W 1 ) compactly. ' 
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Let us show that the operators Qi, ■ ■ ■ ,Q$ are compact. Let v m — > v* weakly in W 2 , T m T * 
weakly in W\. Then v m — > u*, r m — > r* strongly in L q (0,T;W q (ft)) and in L q (0, T; L g (Cl)), and 
Vv m -> Vv*, Vr m -> Vr* strongly in L 9 (0, T; L 9 (fi) n ). 

By Krasnoselskii's theorem [HI [18] on continuity of Nemytskii operators we have 

D { ■ 3 ' 3 f m , T m J > ( ' 3 ' 3 ^* 3 ) 3 

-^1 ' 3 ' ) 3 ) * ^ ( ' 3 '3 ^*3 'T"* ) 3 

strongly in L p ((0,T) x fi) nxn , 

7(-3 -jVm.Tm) -> 7(-, -,V», T*), 

strongly in L p ((0, T) x Q) for all p < 00, and 

f ( ' 3 '3 3 ^"m ) * / ( ' 3 '3 ^*3 ''V ) 3 

strongly in L 2 ((0,T) x fi) n . 

Clearly, if a sequence of functions y m converges in L q ((0,T) x Q), and another sequence z m 
converges in L p ((0,T) x Q), i + ~ = |, then their pointwise products t/ m z m tend to the product of 
their limits in L 2 ((0,T) x Q). 

Hence, D(-,-,v m ,T m )Vv m — > D(-,-,D t) r«)V?J t in L 2 (0, T; L 2 (fi) n ). Therefore 

Qi{v m ,T m ) -> Qi(v*,n) 
in L 2 (0,T; H-\Q)) (and all the more in L 2 (0,T; H- 2 (Q))). Similarly, 

QiiVfni T~rri) * Qi(p* i''*)] ' 2,. ..,5, 

in L 2 (0,T; fl-^fi)). 

Hence, the operator Q is also compact. 
Note that 

(A 2 u,u) = \\u\\l 

for u e H 2 (Q), and 

(A 2 u,u) l = (A 2 u,Au) = (Au,Au) 1 = \\u\\ 2 x 

for u G X . Therefore the operators A\ and A 2 are invertible (e.g. by Theorem 1.1 from [12], Chapter 
VI, or Lemma 3.1.3 from [27]). Hence, A is also (continuously) invertible. 
Rewrite equation (4.26) as 

(u, r) + \A~ x Q{u, r) = A _1 (0, 0, « , So). (4.27) 

A priori bounds from Lemmas 4.1 and 4.2 imply that equation (4.27) has no solutions on the 
boundary of a sufficiently large ball B in W 2 x W\, independent of A. Without loss of generality 
a Q = A (0, 0, u Q , q ) belongs to this ball. Then we can consider the Leray - Schauder degree (see e.g. 
J) of the map / + \A^Q (J is the identity map) on the ball B with respect to the point a , 



deg LS (I + \A-^Q,B, oq). 
By the homotopic invariance property of the degree we have 

deg LS (I + XA~ X Q, B, a ) = deg LS (I, B, a ) = 1 ^ 0. 
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Thus, equation (4.27) (and, therefore, problem (4.1) - (4.4)) has a solution in the ball B for every A. 
□ 

Proof of Theorem 3.1. Take a decreasing sequence of positive numbers e m — > 0. By Lemma 
4.3, there is a pair (v m , r m ) which is a weak solution to problem (4.1)-(4.4) with \ = 1, e = e m . 

Due to a priori estimate (4.6), without loss of generality (passing to a subsequence if necessary) 
one may assume that there exist limits 

u = lim v m , which is *-weak in Loo(0, T; L 2 (f2)) and weak in L 2 (0, T; H 1 (Q)); 

m—+oo 

q = lim r m , which is *-weak in Loo(0, T; H x (f2)) and weak in L 2 (0, T; H 1 (Q)). 

m—roo 

Moreover, due to Lemma 4.2, without loss of generality one may assume that v' m — > u' weakly in 
L 2 (0, T; H^ 2 ), r' m — > q' weakly in L 2 (0, T; H^ 1 ). Then, by [17, Corollary 4], v m — > u, r m — > ? strongly 
in C([0, T]; H^ 1 ). Therefore u and s satisfy (3.8). 

Furthermore, by Corollary 4], v m — > u, r m — > <? strongly in L 2 (0,T; L 2 ). 

By Krasnoselskii's theorem [TH [TS] we have again 

D{-, -,w m ,r m ) -> £>(■, -,u,<;), 
E{-, -,v m ,T m ) -+ E(-, •,«,<?), 

strongly in L p ((0,T) x fi) nxn , 

strongly in L p ((0, T) x f2) for all p < oo, and 

/(-, ;V m ,T m ) -> /(•, •,«,?), 

strongly in L 2 ((0, T) x fi) n . 

Observe that if a sequence of functions y m converges weakly in L 2 ((0,T) x fl), and another 
sequence z m converges strongly in L p ((0,T) x £7), then their pointwise products y m z m converge 

1 + 1 = 1 

P T 2 g- 



weakly to the product of their limits in L 9 ((0,T) x fi), - + | ~ 1 



Therefore, 



D{-, ; V m , T m )Vv m -+D(;;U, $)Vu, 

E(; •, v m , r m )Vr m -»• E(-,-, u, ?)V?, 



weakly in L 9 (0,T; L g (fi) n ), 



7( - ) -i "m, r m )v m -> /?(■, -, u, q)u, 

Pi(; ■,v m ,T m )r m -> ■,«,?)? 

weakly in L q (0,T; L q (Q)) for 1 < g < 2. Therefore the right-hand members of (4.1) converge to the 
corresponding right-hand members of (3.21) weakly in L q (0,T; W~f^(Q)). 

Due to Lemma 4.1, e m \\v m \\ L2 ^,T;H 2 (Q)) = y/^V^\\v m \\L 2 (o,T;H 2 (n)) -> 0. Hence, £ m \\A 2 v m \\ L2 ^ T . H -: 
0. Similarly, £ m ||Tm||L 2 (o,T;X J v) tends to zero, so 
£ m\\A 2 T m \\ L2 ( 0tT , H -i^ — > 0. 

W.l.o.g. we may assume, in addition, that q > Then, by Sobolev theorem, H 2 (Q) C 

Wj/ g _i(^)> so W, -1 ^) C H N 2 (tt). Passing to the limit as m ^ oo in (4.1) and (4.2) with A = 1, 
e = e m , i> = u m , r = r m in the space of distributions on (0,T) with values in H^ 2 (il) (for (4.2) it is 
possible for H^ 1 as well), we conclude that the pair (u, s) is a solution to (3.5)-(3.8). 

It remains to observe that the right-hand side (and, hence, the left-hand side) of (3.21) belongs 
to L 2 (0, T; H^ 1 ), and, due to (3.9), the ones of (3.6) belong to L 2 (0, T; H 1 ). □ 
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5 Long-time behaviour 

Theorem 3.1 implies that solutions to (3.5)-(3.8) can be continued, step by step, onto the whole 
positive semi-axis: 

Corollary 5.1. Given u E L 2 {VL), <Jo £ H 1 ^) and ip G L 2 j oc (0, oo; H]^ (ft)), there is a pair 

u E L 2 , loc (0, oo; H\n)) f| #/ oc (0, oo; H N \n)), q E #,^(0, oo; ^(O)) (5.1) 

which satisfies (3.8), whereas (3.21) holds true in H^^fl) a.e. on (0, oo), and (3.6) holds a.e. in 
(0, oo) x Q. 

Below we keep assuming conditions i)-vi) of Section 3, but we replace (3.17)-(3.19) with stronger 
requirements, namely 
vii) 

\f(t,x,v,r)\ < f(t,x), \g(t,x,v,r)\ < g(t,x) (5.2) 

with some known functions f,gE L 2 ((0, oo) x Q), and 
viii) there ar^j] positive numbers T and Tq such that 



^4 



Z,v) >r (iei 2 + N 2 ) (5.3) 



for any £,,r] E M n . 

Consider any global weak solution (u,q) existing by Corollary 5.1. Denote \& = A~ l ip. Then 
* G L 2t i oc {0, oo; ^(Q)). Assume, in addition, that * G L 2 (0, oo; n Li(0, oo; L 2 (ft)). 

Lemma 5.1. The following estimate is valid: 

IplliooCO.DO^Cfi)) + l|V-u||L 2 (0,oo;L 2 (n)) 
+ I|Vc||l oo (0,oo;Z, 2 (C)) + 1 1 V?| 1 2,3(0,00^3(0)) < C. (5.4) 

Proof. The condition (5.3) can be rewritten as 

p(-)r 2 £,0*» - (j*(-)ti,v)*> + {E(-)r% v ) Rn - (/?(•)£, 77V 

>r (r 2 |£| 2 + M 2 ) (5.5) 

for any 77 G M n (just substitute r£ for £ in (5.3)). 

Take the "bra-ket" of —div^Vqit) E H^^Q) and the terms of (3.6) (as elements of i¥ 1 (f2)) at 
a.a. t E [0,T]: 

(Vc',Vc) = (V(P 1 (t,x,u,q)q + j{t,x,u,q)u),Vq). (5.6) 

Thus, 

^||V?f = {P{t,x,u,q)Vu + ^{t,x,u,q)Vq + g{t,x,u,q),Vq) . (5.7) 

Take the "bra-ket" of (3.21) (as elements of H^ 1 ^)) and u(t) E H 1 ^) at a.a. t E [0, T], arriving 
at (cf. the proof of Lemma 4.1) 

1 d 



7 See [211 Section 5] for a discussion whether this assumption is realistic. 
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Multiply it by T 2 and add this with (5.7): 

= -(£>(*, x, u, ?)r 2 v« + E(t, x, u, <^)r 2 v^, Vu) 

+ x, m, ^) Vw + /i(t, x, m, V?, 

+ x, u, 0, V?) - T 2 (f(t, x, u, 0, Vm) + T 2 (^, «) . (5.9) 
Using (5.5), we conclude that 

< (g(t, x, u, q), V<r) - T 2 (/(t, x, ii, 0, Vm) + T 2 (^, u) i . (5.10) 

Integrating along the interval (0, i), t > 0, we get 

t t 

YKt)ll 2 + ^llv^)|| 2 + r r 2 | ||vM s )|| 2 d s + r 1 ||v ? 00|| 2 <fe 



<?ll«oll 2 + 





t 



2 " u " 2 



^l|V^o|| 2 + y (^(s,a;,u(s),?(s)),V?(s)) 



o 

t 



+ 





I r 2 (*( s )x s ))d s 



+ y [r 2 (v^( s ),v M ( s ))-r 2 (/( s ,x,«( s ),^( s )),v M ( s ))] ds. (5.ii) 



Applying the Cauchy-Buniakowski inequality, Cauchy's inequality and (5.2), we observe that 



t 



J (w(s) + f(s, x, u(s),<;(s)), V«(s)) ds| 



,oo)xQ) + ||/IU 2 ((o,oo)xn)]||V'u|| L2 ((o,t) X Q) 
,oo;L 2 (n)) + ||/||L 2 ((o,oo)xn)] 2 + y ||Vw||| 2(0it . L2(Q)) . (5.12) 

Similarly, 



,oo)xfi) t;i 2 (n))- (5.13) 



And, obviously, 

| y r 2 (M/( S ), W ( S ))rf S | < r 2 ||*|| Ll( o i00;L2( n)) IKII^^W). (5.14) 
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Inequalities (5.11)— (5. 14) yield 

t t 
ylKOI| 2 + ^l|V^)H 2 + ^ J \\Vu(s)\\ 2 ds + ^ J \\V,(s)fds 



< C + C\\u\\ Loo{ o, t .,L 2 (n)), (5.15) 

where C is independent of t, thus we have the same inequality between the essential supremums of 
both members on (0,t), in particular, 

r 2 

y II U \\ L^t-Mm - C + C \\ U \\Loo(0,t;L 2 (n)), 

so 

IMUoo(0,t;L 2 (n)) < C, 

and (5.4) follows from (5.15). □ 

Estimate (5.4) means, in particular, that in a certain sense the concentration u(t) tends to a 
constant as t — > oo. 
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